Topological Defect Density in One-Dimensional 
Priedmann— Robertson— Walker Cosmological Model: Corrections 
Inferred from the Multi-Josephson-Junction-Loop Experiment 



Yurii V. DumiiQ 

IZMIRAN, Russian Academy of Sciences, Troitsk, Moscow region, 142190 Russia 

(Dated: 19 August 2003) 

The data on a strongly-nonequilibrium superconducting phase transition in the multi-Josephson- 
junction loop (MJJL), which represents a close analog of one- dimensional Friedmann-Robertson- 
Walker cosmological model, are used to refine a concentration of topological defects after the phase 
transitions of Higgs fields in the early Universe. The thermal correlations between the phases of order 
parameter revealed in MJJL can reduce considerably the expected number of cosmological defects 
and, thereby, show a new way to resolve the long-standing problem of their excessive concentration. 
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INTRODUCTION 

An important feature of the phase transitions of Higgs 
fields expected at the early stages of cosmological evolu- 
tion is formation of topological defects, whose concentra- 
tion can be roughly estimated as n ~ 1/ Ccff • Here, d — i, 
2, and 1 for the monopoles, cosmic strings, and domain 
walls, respectively; and S^cS is the effective correlation 
length, which is commonly assumed to be less than the 
cosmological causality horizon, ^off ^ c / H^t , where c is 
the speed of light, and iJpt is Hubble constant at the in- 
stant of phase transition. So, the resulting concentration 
of the defects should be 

n>{H^tlc)\ (1) 

which exceeds considerably the upper limits following 
from the observational data (e.g. review 

For the sake of definiteness, we shall further consider 
the simplest type of defects — the domain walls. In such 
a case, their excessive concentration can result in the t^- 
dependence for the evolution of cosmological scale factor, 
leaving less time for galaxy formation and changing the 
rate of nucleosynthesis. In addition, the domain walls 
can produce unreasonably large distortions in the cosmic 
microwave background radiation [3, 01 ■ 

A commonly-used approach to resolve the domain-wall 
problem is to introduce some mechanism of their annihi- 
lation, for example, utilizing the concept of the so-called 
"biased" (or asymmetric) vacuum |2|. As a result, under 
appropriate choice of parameters of the fields involved, 
the regions of "false" (energetically unfavorable) vacuum 
will quickly disappear, eliminating the corresponding do- 
main walls. (A detailed discussion of the various regimes 
of evolution can be found in Ref. 01-) 

Unfortunately, the concept of the biased vacuum is 
not sufficiently supported by realistic models of elemen- 
tary particles. So, it becomes interesting to seek for so- 
lution of the problem by another way, namely, to an- 
swer the question if there are some mechanisms reduc- 
ing the commonly-used lower limit on the topological 



defect concentration To examine the efficiency of 

formation of the defects by the strongly-nonequilibrium 
symmetry-breaking phase transitions, a number of ex- 
periments was carried out recently in the condensed- 
matter systems, such as liquid crystals |^|^, supcrfluid 
''He m and '^He 0, S, ^'^'^ superconductors both in the 
form of bulk samples |9| and quasi-one-dimensional struc- 
tures [loj . (The comprehensive review was given, for ex- 
ample, by Kibble ^.) 

As was established in the last- mentioned experiment, 
utilizing the multi-Josephson-junction loop (MJJL) ^lo| . 
the probability of occurrence of various field configura- 
tions during the phase transition can be estimated much 
better if the energy concentrated in the defects is taken 
into account. The main aim of the present article is to 
apply this idea to the consideration of cosmological phase 
transitions and, thereby, to find the range of parameters 
at which formation of the domain walls can be substan- 
tially suppressed even without any assumption of asym- 
metry (or bias) of the vacuum states. 

REVIEW OF THE MJJL EXPERIMENT 

A general scheme of the MJJL experiment is shown in 
Fig. ^ a thin loop produced of the high-temperature 
YBa2Cu307_5 superconductor and interrupted by 214 
Josephson junctions at the grain boundaries is rapidly 
cooled from the normal to superconducting phase 
(namely, from -100 K to 77 K) 

At Tc = 90 K, the segments separated by junctions 
become superconducting; however the junctions are still 
normal and, therefore, the superconducting segments are 
effectively separate. So, a random phase of the supercon- 
ducting order parameter should be established in each of 
them. 

At subsequent cooling down to the temperature Tcj 
(which is 5-^7 K below Tc), the Josephson junctions also 
become superconducting, so that a superconducting cur- 
rent / is formed along the entire loop. As a result, the 
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FIG. 1: Sketch of the MJJL experiment. 



loop will be penetrated by the magnetic flux which is 
just the measurable quantity. (This is quite close to the 
original idea by Zurek "^^o proposed to observe a 
spontaneous rotation produced by the superfluid phase 
transition in a thin annular tube.) 

Because of the phase jumps between the isolated seg- 
ments formed at the stage when T^j < T < Tc, the final 
magnetic flux through the loop is nonzero and varies ran- 
domly from one heating-cooling cycle to another. Its his- 
togram A'$($) derived from 166 cycles is well described 
by the normal (Gaussian) law with zero average value 
and standard deviation 7.4 00 (where (p^ is the magnetic 
flux quantum) . 

In fact, the above- written experimental value is unrea- 
sonably large: if the phase jumps between the segments 
would be absolutely independent of each other, then the 
expected width of the distribution should be only 3.6 0o- 
On the other hand, the excessive experimental value was 
satisfactorily explained by the authors of this experiment 
under assumption that the phases of superconducting or- 
der parameter in the isolated segments are correlated to 
each other so that probability P{Si) of phase jump Si at 
the i'th junction is given by Gibbs law: 



P{S,) cx eM-EJiS^)/kBT] 



(2) 



where Ej is the energy concentrated in the Josephson 
junction, T is the temperature, and is Boltzmann 



constant. 

So, the main conclusion from the results of the above 
experiment is that the energy concentrated in the defects 
should be taken into account in calculating the proba- 
bility of realization of the corresponding field configura- 
tions, even if the defects are located at distances exceed- 
ing the effective correlation length of the phase transi- 
tion. In other words, the correlation length represents a 
size of the minimal domain in which the field must be 
uniform, but it cannot be considered as the scale beyond 
which the field states are absolutely independent of each 
other. Of course, for the correlations at the larger scales 
to take place, the corresponding regions should have the 
possibility to interact with each other in the course of 
the previous evolution, before the phase transition. (For 
the sake of brevity, we shall call them to be in a coherent 
state.) The last- mentioned condition is satisfied auto- 
matically in the laboratory systems but requires a spe- 
cial consideration in the cosmological applications (see 
inequality H1U|I below). 



COSMOLOGICAL IMPLICATIONS 

First of all, we should emphasize close similarity be- 
tween the symmetry-breaking phase transitions in the 
multi-Josephson-junction loop, drawn in Fig. ^ and the 
simplest one-dimensional (ID) Friedmann-Robertson- 
Walker (FRW) cosmological model. To make the quan- 
titative estimates, let us consider the space-time metric 



ds'^ = dt'^ - a^{t) dx^ 



(3) 



(where a{t) is the scale factor of FRW model) and the real 
scalar field ip (simulating Higgs field) whose Lagrangian 

L{x,t) = i [{dt^f- (9.^)2] - ^ [vp^- (^VA)]' (4) 

possesses Z2 symmetry group, to be broken by the phase 
transition. (The resulting topological defects in such a 
model will be evidently the pointlike domain walls.) 

As is known, the stable low-temperature vacuum states 
of the field are 



(5) 



while a domain wall between them is described as Lp [x) — 
± tanh[(/i/\/2 ) [x — xo)] and involves the energy 



E ■ 



2V2 fi^ 



(6) 



(From here on, it will be assumed that thickness of the 
wall is small in comparison with a characteristic 

domain size.) 
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FIG. 2: Conformal diagram of the space-time. 



Next, by introducing the conformal time rj — J dt/a{t) , 
the space-time metric (0) can be reduced to the confor- 
mally flat form p^ : 



(7) 



so that the Ught rays {ds^ = 0) will be described by the 
straight lines inclined at ±7r/4 : x = ztrj + const. 

The entire structure of the space-time can be conve- 
niently presented by the conformal diagram in Fig. |21 
Let 77 = and ^7 = f?o be the beginning and end of the 
phase transition, respectively, and rj — r]^ be the instant 
of observation. Then, as follows from a simple geometric 
consideration, 770 is the maximum correlation length; so 
that 

N ^ {ri^~ rio)/m ^ V*/Vo (at large iV) (8) 

is the minimum number of spatial subregions causally- 
disconnected during the phase transition. (Their final 
vacuum states are arbitrarily marked by the arrows.) 

A probability of phase transition without formation of 
the domain walls is commonly estimated as ratio of 
the number of field configurations without domain walls 
to their total number: 



Pn = 2/2 



N 



(9) 



which tends to zero very sharply at large N. So, the ob- 
servable part of the Universe, represented by the large 
upper triangle in Fig.|21 will inevitably contain a consid- 
erable number of the domain walls. 

On the other hand, if a sufficiently long interval of the 
conformal time 



A77 > 77* 



(10) 



preceded the phase transition, then a coherent state of 
the Higgs field (shown by the lower shaded triangle) will 



be formed by the instant 77 = in the entire region ob- 
servable at ?7 = r/^ . 

The inequality 1)1 0|l can be satisfied, particularly, in 
the case of sufficiently long de Sitter stage. Really, if 
a{t) = exp{Ht), then 



1 -Hf 

77 = — — e + const 
H 



—00 at t 



(11) 



so that A77 can be sufficiently large. From the viewpoint 
of elementary-particle physics, the de Sitter stage can be 
easily realized in the overcooled state of Higgs field just 
before its first-order phase transition. (Let us remind 
that just this idea was the basis of the first infiationary 
models; for more details, see review 0|.) 

Next, if condition IjlOII is satisfied, it is reasonable to 
assume that the coherent state of the field (p will exhibit 
Gibbs-like correlations (similar to the ones occurring in a 
superconducting Bose condensate of MJJL) between the 
all N subregions drawn in Fig. [3 In such a case, the 
probability should be calculated taking into account 
the Gibbs factors for the field configurations involving 
domain walls: 



pO 



2/Z, 



(12) 



E E exp{-f (13) 



where 



=±1 

Here, Sj is the spin-like variable describing a sign of the 
vacuum state in the j'th subregion, E is the domain wall 
energy, given by ®, and T is the characteristic temper- 
ature of the phase transition. 

From the mathematical point of view, statistical 
sum (|13|) is very similar to the one for Ising model, well 
studied in the condensed-matter physics (e.g. Ref. [l5|). 
Using exactly the same methods (attention should be 
paid to the appropriate choice of zero energy), we get 
the final result: 



pO 



[1 



-E/T' 



N 



[1 



,-EIT' 



N 



(14) 



(Yet another method of calculating this quantity, based 
on explicit expressions for the probabilities of field config- 
urations with various numbers of the domain walls, was 
presented in our article 

As is seen in Fig. |31 P^ drops very sharply with in- 
creasing N at small values oi E/T |21| but becomes a 
gently decreasing function of N when the parameter E/T 
is sufficiently large. Therefore, just the large energy con- 
centrated in the domain walls turns out to be the factor 
substantially suppressing the probability of their forma- 
tion. 

As can be easily derived from H14|l . the probability of 
absence of the domain walls in the observable region of 
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FIG. 3: The probability of phase transition without forma- 
tion of the domain walls as function of the number of 
disconnected subregions A'^ and the ratio of the domain wall 
energy to the phase transition temperature E/T. 



space-time becomes on the order of unity (for example, 
1/2) if E/T > In N or, by substituting © and ©, 



> In^, 
AT 770 



(15) 



which is just the required inequality, mentioned in Intro- 
duction. Because of the very weak logarithmic depen- 
dence in the right-hand side of 1)1511 . this condition may 
be quite reasonable for a certain class of field theories. 

Moreover, the situation should be even more favorable 
in the cases of higher spatial dimensions. A well-known 
property of Ising models in 2 and 3 dimensions is the ten- 
dency for aggregation of the domains with the same sign 
of the order parameter when the temperature drops be- 
low some critical value (e.g. Ref. 0|). In the condensed- 
matter applications, this corresponds to the spontaneous 
magnetization of a solid body. Regarding the cosmologi- 
cal context, we can expect that probability of formation 
of the domain walls will be reduced dramatically at the 
sufficiently large values oi E/T 22]. This sub ject will be 
discussed in more detail in a separate paper |l8|| . 

Although the model under consideration requires a 
presence of inflationary stag e (for the coherent state of 
Higgs field to be formed 23]), it is not equivalent to the 
traditional inflationary models, where decrease in con- 
centration of the previously formed topological defects is 
a purely geometrical effect, associated with sharp expan- 
sion of the space. This raises a problem of new topolog- 
ical defects that can be formed at a subsequent stage of 
decay of the inflaton field by itself. On the other hand, 
the inflationary stage in our model results in suppress- 
ing formation of the defects (rather than subsequent de- 
crease in their concentration). From this point of view, 
it represents a self-consistent solution of the problem. 
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The actual loop used in the experiment is a winding strip 
engraved at the boundary between grains of the super- 
conductor film. 

The standard deviation of the total distribution is just 
a typical value of the magnetic flux measured in each 
individual heating-cooling cycle. 

E/T = corresponds to the case when Gibbs-like corre- 
lations between the subregions are absent at all, so that 
formula (1141 coincides with @. 

To avoid misunderstanding, it should be emphasized once 
again that the Ising model considered here is only an 
auxiliary mathematical construction, describing a final 
distribution of the domain walls after the phase transi- 
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tion. This physical phase transition of the field bears easily if the phase transition develops from a strongly 

no relation to the phase transition in Ising model. overcooled state of the Higgs field. 

[23] Besides, the inequality (1151 can be satisfied much more 



